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We report on numerical results from a revised hydrody- 
namic simulation of binary neutron-star orbits near merger. 
We find that the correction recently identified by Flanagan 
significantly reduces but does not eliminate the neutron-star 
compression effect. Although results of the revised simula- 
tions show that the compression is reduced for a given total 
orbital angular momentum, the inner most stable circular or- 
bit moves to closer separation distances. At these closer or- 
bits significant compression and even collapse is still possible 
prior to merger for a sufficiently soft EOS. The reduced com- 
pression in the corrected simulation is consistent with other 
recent studies of rigid irrotational binaries in quasiequilibrium 
in which the compression effect is observed to be small. An- 
other significant effect of this correction is that the derived 
binary orbital frequencies are now in closer agreement with 
post-Newtonian expectations. 

PACS number(s): 97.80. Fk, 04.25.Dm, 04.40.Dg, 97.60. Jd 



I. INTRODUCTION 

In a recent paper E. Flanagan JlJ] has pointed out an 
inconsistency in the solution of the shift vector in our 
previously reported numerical hydrodynamic simulations 
of binary neutron stars in quasiequilibrium orbits. 
In that paper it was suggested that this may be the source 
of the controversial relativistic hydrodynamic result that 
as the stars approach each other their interior density is 
observed to increase during numerical simulations. 

In view of the controversial nature of the compression 
effect H and the fact that no other relativistic hydrody- 
namic treatments are yet capable of attacking this prob- 
lem, it is of course essential to incorporate this correction 
into the previous hydrodynamic simulations to see if in- 
deed this is the source of the observed compression effect. 

In this short note we report of the first results of sim- 
ulations in which this correction has been applied. We 
find that the correction does not eliminate the compres- 
sion effect. It does, however, significantly diminish its 
magnitude at a given angular momentum. Perhaps more 
importantly, the correction causes the orbital frequency 
of the stars to increase such that the orbits are in closer 
agreement with post-Newtonian expectations. At the 
same time, however, this causes the inner most stable 
circular orbit (ISCO) to move to closer separation dis- 



tances and higher associated fluid velocities. The com- 
pression continues to scale with the magnitude U of the 
spatial component four velocity as noted in Jj^]. Hence, 
by the time of inspiral and merger a significant com- 
pression effect is still possible. For a sufficiently soft 
equation of state, neutron stars in the mass range of ob- 
served neutron-star binaries still collapse to individual 
black holes. 

II. THE CORRECTION 

The correction pointed out by Flanagan concerns the 
the solution for the space-time components of the ADM 
(3+1) metric 

ds 2 = -(a 2 - (3i(3 l )dt 2 + 2(3idx l dt + j.jdx'dx 3 , (1) 

where the space-time components of the metric, ft 1 are 
referred to as the shift vector. 

We solve for the components of the shift vector by 
applying the ADM momentum constraint Q] , 

Di(K ij -<y ij K) = 8nS J . (2) 

Where Dj is the three-space covariant derivative |Q , and 
S l is the ADM three-momentum density. 

In our approximation scheme, we impose maximal- 
slicing [Tr(iT y ) = 0] and demand that the spatial three 
metric 7y be conformally-flat. Under these constraints, 
the second term on the left hand side of Eq. vanishes 
and we have, 

D t K 13 = 8TrS j . (3) 

Ultimately, we can reduce equation (^) to a Poisson-like 
equation for the shift vector, 

W =4^-^(^-/3) . (4) 

The correction pointed out in [jlj is in the derivation of 
the source density pp. In la] the contravariant spatial 
components of the four-momentum density were utilized 
in place of the contravariant three-momentum density. 
[This is valid for the covariant components but not for 
the contravariant components.] Hence, in J3j the source 
density was incorrectly written as, 
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1 d£, fdp dp 2 dp\ 



(5) 



where a = p + pe + P is the inertial mass density, £ 
a/</> 6 , = at/*, and 

Si=aWUi . 



(6) 



The correct source should not have the second term 
on the right hand side. This is a spurious term which 
came from an erroneous identification of the contravari- 
ant momentum density with the ADM contravariant 
three-momentum density. The correct source density is, 

* ^4C, 1 ^ (dp* dp 2 d{3 k \ 

where Si is still given by Eq. (Q). 

This correction would have little effect on the final re- 
sults if the spurious 4pW 2 a term were small compared 
to the term containing $j. That is not the case how- 
ever. In the simulations of ||] it was observed that those 
two terms nearly cancelled leading to small values for the 
p. With this cancellation no longer in effect, the source 
for the shift vector equation is now a substantial quan- 
tity. This has several implications for the orbit dynamics 
which we now analyze. 



III. RESULTS 

A. r = 2 EOS 

As one way to identify the effect of this correction we 
have done calculations similar to the bench-mark calcu- 
lation of ||, i.e. we have employed a simplistic T = 2 
polytropic equation of state (EOS), P = Kp r , where 
K = 1.8 x 10 5 erg cm 3 g~ 2 . This gives a maximum 
neutron-star mass of 1.82 Mq. The gravitational mass 
of a single tub — 1.625 M Q star in isolation is 1.51 Mq 
and the central density is p c = 5.90 x 10 14 g cm -3 . The 
compaction ratio for an isolated star with the same nu- 
merical grid resolution (in terms of Schwarzschild coor- 
dinates) is m/R = 0.14. [Note that the central density 
for these stars these stars is slightly higher than quoted 
in H (for which p c = 5.84 x 10 14 g cm -3 ). This is due 
to slight changes in the finite differencing of the present 
calculation.] The binary is taken to have a fixed angu- 
lar momentum J = 2.5 x 10 11 cm 2 (J/Mq — 1.27 where 
Mg = 2mc). In quasicquilibrium circular orbit for this 
angular momentum, the proper (coordinate) separation 
distance between centers is 118 (102) km with a frequency 
consistent with post-Newtonian estimates. The stars re- 
lax to a nearly irrotational flow, and the central density 
only slightly increases by 0.8% to 5.95 x 10 14 g cm~ 3 as 
summarized in Tabic [|. 



This very slight increase in central density is compara- 
ble to the numerical accuracy of the hydrodynamic cal- 
culation which we estimate to be ~ ±0.5%. It is also 
consistent with the very small amount of compression 
noted in simulations of purely irrotational T — 2 stars 



by several groups 
tio. For example, in 



with a similar compaction ra- 
stars with a similar compaction 
ratio exhibit an increase in central density of 0.1% at 
this separation distance. In |sj changes in central den- 
sity of order 0.1 ± 0.5% are consistent with the numerical 
results; and in |10|, stars with m/R — 0.14 show a max- 
imum central compression of about 2% at a coordinate 
separation distance of cI/Rq = 1.6 in their notation. A 
reasonable extrapolation of their figure 6 to our separa- 
tion {d/R w 4.5) suggests an increased central density 
consistent with our results. 



B. Realistic EOS 

In all of the above simulations [|]-[l0) ■ it has been noted 
that the compression effect dramatically increases for an 
EOS with an increased compaction ratio and closer or- 
bits. Hence, in the remaining discussion we consider the 
"realistic" neutron star equation of state from for 
which the compaction ratio is much higher m/R « 0.25. 
For this set of calculations, we not only compare with 
previous published results , but also quantify the ef- 
fect of the correction on the location of the ISCO and the 
orbital frequency for the binary. 

Consider first the somewhat soft EOS with a critical 
mass of m c = 1.575 and stars with a baryon mass of 1.548 
M corresponding to M G = 1.39 M p c = 1.34 x 10 15 g 
cm -3 in isolation. These parameters roughly correspond 
lO to one of the simulations in Table 2 of Ref. {|. Table 
y compares differences between the previously derived 
quantities B and the corrected values for the J = 2.7 x 
10 11 cm 2 run. For the uncorrected runs, this J was the 
last orbit before the stars collapsed. 

Here, similar trends to those noted in the r = 2 EOS 
are to be noted. The correction significantly stabilizes 
the stars for this J. However, the amount of compression 
(~ 3%) for this EOS is greater than for the simulation in 
Table § 



In Table 



III we show properties of corrected orbits as 



the stars approach the ISCO. For the uncorrected sim- 
ulations the ISCO for these stars with a stiff EOS (i.e 
no collapse) occurred at J « 2.0 x 10 11 cm 2 . In the 
corrected runs the orbital frequencies are larger for the 
same fixed angular momentum. This moves the ISCO 
in to closer distances and smaller J. Thus, although the 
compression effect is less for a given J, the stars reach 
higher velocities before the ISCO. This allows for signif- 
icant compression (~10% in this case) before the ISCO. 
Note, however, that these stars do not collapse as they 
had done in the uncorrected simulations. 
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We did find, however, that collapse could occur if the 
stars were increased in mass from m,Q = 1.39 to 1.44 
M Q (rriB = 1.61 M Q ) for this simulation. The results 
from this run are given in Table IV . Collapse of the stars 
was observed to occur for very close orbits (J = 1.85 x 
10 11 cm 2 ) just before inspiral. The coordinate separation 
between stars was only 2.4 times the coordinate radii. At 
,/ = 1.8 x 10 11 cm 2 the orbit still appeared stable though 
the stars wanted to collapse. In another calculation we 
used a softer EOS for which m c = 1.54 M Q and ttiq — 
1.40, mb = 1.54 M Q . For this case, collapse occurred 
with J = 2.0 x 10 11 cm 2 . 

Thus, collapse may still be a possibility albeit for stars 
close to the maximum mass of a soft EOS and for very 
close orbits. Such is soft EOS is a reasonable possibility. 
For example, collapse would always occur prior to inspiral 
for typical-mass neutron stars modeled with the EOS of 
Bcthc and Brown (l2| . 



velocity. This remains true even though the magnitude 
of U for the binary slso slightly decreases. 

A way of understanding the decrease in U is to con- 
sider the balance between the relativistic analog of the 
centrifugal force and the gravitational force M. The j3 
part of the centrifugal force can be written as Sjdftf /dx % . 
This centrifugal force consists of two parts. One part 
scales as U 2 and one scales as (3U. Since the (3 term is 
now much larger the magnitude of U required to balance 
the gravitational force becomes smaller. 

Alternatively the reason for a decrease in U while uj 
increases can be seen from the definition of angular mo- 
mentum. In our simulations we constrain each orbit to a 
specific value for the covariant z component of the orbital 
angular momentum vector. In the corrected simulations 
the shift vector contributes significantly to the covariant 
angular momentum, 



J = J; 



J 



(9) 



C. Angular Momentum at the ISCO 

Another significant effect concerns the specific angular 
momentum J/M G as the stars approach the final orbits. 
In the previous results |l3| the specific angular momen- 
tum at the ISCO was sig nificantly higher (J/M G « 1.3) 
than that of a maximally rotating Kerr black hole. This 
would imply complicated dynamics during inspiral be- 
fore the stars could merge. For the corrected results, 
J/Mq — 1.03 for the last computed stable orbit of the 
sequence in Table |l| For J/M G = 0.99 the orbit is 
unstable. Thus, we expect that when the stars begin 
to inspiral the specific angular momentum is very near 
unity, J/Mq w 1.00 and will become < 1 as the orbit 
plunges. Hence, the stars can immediately spiral inward 
to form a Kerr black hole near maximum rotation. This 
has important implications for the emergent gravity wave 
signal from the subsequent ringing. 



IV. DISCUSSION 

Perhaps the most significant effect concerns the orbital 
frequency. The correction has caused the frequency to in- 
crease to be much closer to the expected post-Newtonian 
result (cf. Table |J). It is easy to understand why this 
is so. In the hydrodynamic calculations, the orbital fre- 
quency uj is determined by minimizing the average coor- 
dinate three velocity, i.e. 



(V) oc (U - (3 - uj x r) « 



(8) 



The vector quantity is of opposite sign to U and is 
much larger in the corrected simulations. This means 
that a larger value for ui is required to minimize the three 



where 



and 



where 



yS x — d x , 
a 



J? 



0/3 



y(3 x - x(3 v 



0/3 



W 2 a -pe-P 



W a 



(10) 



(11) 



(12) 



In the corrected equations Jg is a significant fraction of 
the total angular momentum (<l/2, see Table ||). Thus, 
the corresponding momentum density contributions (J^) 
for a given total angular momentum is less. This con- 
strains U to be smaller in equation 

In it was noted that the increase in central density 
scales with the four velocity for the binary as p c oc U 4 . It 
is note worthy that the U 4 scaling with a similar propor- 
tionality constant still applies to the corrected results as 
shown in Table [V. The main difference, however, is that 



the values of U are lower after the correction. Hence, the 
compression is less for a given total J. 



V. CONCLUSION 

In summary, we have applied the revision recently 
pointed out by E. Flanagan |lj into a hydrodynamic sim- 
ulation of binary neutron stars near the ISCO. The ef- 
fects of this correction have been analyzed. Although 
the magnitude of the compression effect is reduced for a 
given fixed total angular momentum, some compression 
effect remains. The scaling of the compression with the 
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spatial components of the four velocity U remains the 
same. The main difference is that U is smaller for a fixed 
angular momentum. This is at least in part because a 
significant fraction of the covariant angular momentum 
now arises from the corrected shift vector term. 

Another outcome is that the corrected orbital angu- 
lar frequencies are higher than the uncorrected results. 
This brings the anticipated orbital and gravity wave fre- 
quences into closer alignment with post-Newtonian es- 
timates. Another consequence of the higher orbital an- 
gular frequencies is that the ISCO for the binary moves 
in to closer separation distances, higher velocities, and 
higher orbital angular momenta than previously esti- 
mated. This increases the compression such that col- 
lapse to two black holes is still possible before inspiral, 
although only for a softer EOS and/or higher mass. The 
possibility that the EOS is this soft is not yet ruled out 

Hug. 

Although we note here that some compression remains 
in the revised hydrodynamics, its magnitude is consider- 
ably reduced. The question remains then as to whether 
the remaining effect is real or an artifact of the uncertain- 
ties HD introduced by the conformally flat condition on 
the metric. The resolution of this question, however, will 
require simulations in which the full Einstein dynamics 
are included in the orbits. We are currently developing 
a perturbation expansion to examine this question. 
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TABLE I. Comparison between corrected and uncorrected 
results for m_g = 1.625 M (mc = 1.51 M in isolation) stars 
in a binary with J = 2.5 x 10 11 cm 2 and a T = 2 EOS. 



Quantity 


Isolated 


Corrected 


Uncorrected 




Star 


Binary 


Binary 


dp (km) 


oo 


118 


138 


p c (10 14 g cm" 3 ) 


5.90 


5.95 


6.68 


(Sp/po) 


0. 


0.008 


0.14 



TABLE II. Comparison between corrected and uncor- 
rected results (from [4]) for m B = 1.548 M (tug = 1-39 
Mq in isolation) stars in a binary with J = 2.7 x 10 11 cm 2 
and a realistic EOS. 

Quantity Corrected Uncorrected 



TABLE IV. Summary of corrected results for m B = 1.61 
M (m,G = 1.44 M in isolation) binary stars with a realistic 
EOS. 

J (10 11 cm 2 ) lo (rad sec -1 ) dp (km) p c (10 15 g cm" 3 ) U 2 



oo 


0.0 


oo 


1.38 


0.0 


3.0 


1710 


56.4 


1.51 


.013 


2.6 


1800 


52.0 


1.55 


.015 


2.4 


2100 


47.2 


1.57 


.018 


2.2 


2400 


44.2 


1.60 


.020 


2.0 


3000 


37.2 


1.64 


.030 


1.9 


3600 


29.8 


1.72 


.040 


(Collapsing) 


4500 


23.6 


4.05 


.070 



dp (km) 
p c (10 15 g cm" 3 ) 


61 


68 


1.38 


1.98 


(Sp/po) 


0.03 


0.51 


lo (rad sec- 1 ) 


1330 


732 


(uj/lopn) 


1.01 


0.66 


U 2 = W 2 - 1 


0.0091 


0.034 


Jl (10 11 cm 2 ) 


1.4 


2.7 


JS (10 11 ) cm 2 


1.3 






0.55 


.44 



TABLE III. Summary of corrected results for mj = 1.548 
Mq (mc = 1.39 Mq in isolation) binary stars with a realistic 
EOS. 

J (10 11 cm 2 ) lo (rad sec -1 ) d p (km) p c (10 15 g cm" 3 ) U 2 



oo 


0.0 


oo 


1.34 


0.0 


3.0 


1025 


72.0 


1.37 


.0080 


2.8 


1190 


65.6 


1.375 


.0085 


2.6 


1460 


56.2 


1.38 


.0097 


2.4 


1660 


51.4 


1.40 


.0115 


2.2 


2110 


43.8 


1.43 


.0145 


2.1 


2425 


42.0 


1.44 


.0165 


2.0 


2530 


41.2 


1.45 


.0185 


1.9 


2750 


39.0 


1.455 


.020 


1.8 


3000 


35.0 


1.46 


.023 


1.7 


3450 


34.0 


1.47 


.0245 


1.65 


4200 


29.2 


1.50 


.0335 


1.6 


(Inspiral) 
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